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Using concepts from fluid dynamics, a partial differential equation for the shape evolution of a 
metallic nanowire is derived from a semiclassical energy functional that includes electron-shell 
effects. A rich dynamics, involving movement and interaction of kinks connecting locally stable 
radii, leads to the formation of a wire whose equilibrium shape is universal and consists of a 
cylindrical part connected to unduloid-like leads. The universality of the equilibrium shape 
may provide an explanation for the formation of cylindrical nanowires observed in recent 
experiments. 



1. INTRODUCTION 

Recent transmission electron microscopy (TEM) experiments by the groups of Takayanagi^"^ 
and Ugarte^'^ have found that long, almost perfectly cylindrical Au and Ag nanowires are 
frequently formed and have lifetimes of the order of seconds at room temperature. In those 
experiments, holes are burned through a thin metallic film using high-intensity electron beams 
until two holes come close together, leaving a thin, initially rough nanowire in-between. After 
reduction of the electron irradiation to turn into imaging mode, the wires are found to evolve 
towards almost perfect, long cylinders. Subsequent thinning of the wire is found to occur 
through motion of a "kink" from one end of the wire to the other. ^ 

The stability of long cylindrical wires might seem puzzling since nanowires have a large 
surface to volume ratio, suggesting that a wire longer than its circumference should break 
up under surface tension (Rayleigh instability^). This apparent paradox is resolved theoreti- 
cally^"^° with the inclusion of quantum effects, which have been shown to stabilize the wires 
up to large temperatures for a set of "magic" radii. The stability arises through a competition 
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of shell-effects, comparable to what happens in metal clusters, and an interplay of Rayleigh 
and Peierls instabilities. ^° 

A dynamical model, ^^ including surface self-diffusion of a continuous ionic background, 
indicates that cylindrical wires form spontaneously from random initial wires. In this paper, 
the dynamical model is described and applied to the evolution of initially random nanowires, 
emphasizing striking similarities between simulations and available experimental results. 

2. MODEL 

We use a model of a nanowire consisting of free electrons confined to an axisymmetric 
wire by a hard- wall potential. ^^ A wire parallel to the z-axis is characterized by its radius 
R{z,t) in the interval [0,L], and periodic boundary conditions'^ are used to extend R{z). 
The restriction to axisymmetric wires, though it ignores the stability of a few Jahn- Teller 
distorted wires at low conductance,'^ is not expected to alter significantly the dynamics of 
nanowires, since the most stable nanowires have been shown to be cylinders, and axisymmetry, 
once present, is preserved by the diffusion equation, so that Jahn- Teller distortions are most 
likely suppressed dynamically. The free-electron model requires good screening and a spherical 
Fermi surface, both conditions which are met by alkali metals, such as sodium, and to some 
extent by noble metals. Throughout this paper, we use parameters corresponding to Na, but 
results should be qualitatively unchanged for Au and Ag, which are more frequently used in 
imaging experiments. 

In the spirit of the Born-Oppenheimer approximation, the total energy of the nanowire is 
taken to be the electronic energy. Since we are dealing with an open system of electrons, the 
grand-canonical potential is used, and can be separated into Weyl and mesoscopic contribu- 
tions.'^''^ Dropping the volume contribution, assumed to be constant so that the volume per 
atom is fixed, the energy functional is 

n[T,R{z,t)]=aiT)S[Riz,t)]+ [ V{T, R{z,t))dz, (1) 

Jo 

where a = ^^^ is the surface tension,'^ ep and kp are respectively the Fermi energy and 
wavevector, S is the surface area of the wire, and l^ is a mesoscopic electron-shell potential. 



Higher-order terms^^'^^ proportional to the mean curvature, etc., are neglected. The electron- 
shell potential V can be expressed in terms of a Gutzwiller-type trace formula^^ 
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where the sum includes all classical periodic orbits (f , w) in a disk billiard, ^^ characterized by 
their number of vertices v and winding number w, L^^ = 2vRsm{iTw/v) is the length of orbit 
{v,w), and 9^yu = kpL^w — 3f7r/2. The factor f^w = 1 for t> = 2w, f^w = 2 otherwise, accounts 
for the invariance under time-reversal symmetry of some orbits, and a^wiT) = t^w/ ^i^^Tvw 
{tvw = T^kpL^yjT /2Tf-, Tp being the Fermi temperature) is a temperature-dependent damping 
factor. A semiclassical formula similar to Eq. (2) was obtained for perturbations of a cylinder, ^^ 
but we point out that Eq. (2) remains valid for large deformations, as long as new classes of 
non-planar orbits can be neglected (adiabatic approximation).^^ 
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Figure 1. Mesoscopic electron-shell potential V{R,T), for three different temperatures. 



Eqs. (1) and (2) define an energy functional which is general and simple enough to solve for 
nontrivial nanowire geometries with a wide range of radii. The multiple deep minima of V{R), 
shown in Fig. 1, favor certain magic radii^^ and suggest a multistable field theory analogous to 



the sine-Gordon model. Solutions consisting of cylindrical segments connected by kinks may 
therefore be anticipated. 

The dynamical model is derived from two main assumptions: First, in the thin wires con- 
sidered, a significant fraction of the atoms are at the surface and thus surface self-diffusion is 
the dominant mechanism of ionic motion. Second, under the Born-Oppenheimer approxima- 
tion, the electronic energy (1) acts as a potential for the ions, and the dynamics derive from 
ionic mass conservation: 

dR{z,t) Va 9 r . 

where Va = Sn'^/kp is the volume of an atom, and the 2;-component of the surface current is 

given by Pick's law: 

^ PsDs 1 _ g/i . . 

where R' = dR/dz, and ps and Ds are the surface density of ions and the surface self- 
diffusion coefficient, respectively. The precise value of Ds for alkali metals is not known, but 
it can be removed from the evolution equation by rescaling time to the dimensionless variable 
r = {psDsTF/T)t. The chemical potential p, calculated as the change in the energy (1) with 
the addition of an atom at point z, is 
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where KwiT) = a^w coshr^^^,. Eq. (3) is solved numerically using an implicit scheme. 

This model, though neglecting the atomic structure of the wires, is expected to give rea- 
sonable insights into the dynamics and equilibrium shape of metallic nanowires: Delocalized 
electrons determine the cohesion of metals, so that the confinement of the electrons to the wire 
should affect the geometry of the wire in a significant manner. To lowest order, atoms may 
be expected to arrange themselves in order to fit in the shape determined by the electrons. 
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Figure 2. Evolution of an initially (top) random wire of length kpL = 600. Bottom wire 
shows the equilibrium shape, while the middle wire shows an intermediate stage of evolution. 
The conductance and dimensionless time of evolution are indicated above each wire. 



3. EVOLUTION OF A RANDOM WIRE 

In order to mimic TEM experiments/^^ we perform simulations starting from an initially 
random wire given by 



i?(^,t = o)= J2 ^9e''?^ 



(6) 

where q = 27m/L with n an integer, Qmax is of the order of the Fermi wavevector kp, and 
R* = R-q to ensure the reality of R{z). The coefficients Rq and R±i are chosen so as to set 
the initial average radius and curvature of the wire, while other coefficients are random. We 
performed 14 simulations, starting with different random wires of various radii and lengths. 
Figure 2 shows such a wire, whose evolution is typical, at three stages of evolution: (i) Initial 
random shape; (ii) After a relatively short time, the short wavelength roughness has been 
smoothed out, leaving a few cylindrical segments, connected by kinks; (iii) Eventually, all 
kinks move towards the "leads", yielding an equilibrium shape consisting of a cylindrical part, 
whose radius and length depend on the initial conditions, and thicker contacts, whose shape 
will be discussed in the next section. This suggests that the natural evolution of a nanowire. 



at a temperature sufficient for atoms to diffuse, is to form a cylinder, providing an explanation 
of the TEM observation. 

Figure 3 shows details of the thinning process: The top solid line shows the central part 
of a wire, with a conductance G = 8Gq, right before the thinning starts with the nucleation 
of a kink/antikink pair close to one end of the wire (dashed line) . One of the kinks is quickly 
absorbed into the closest lead, while the other moves at constant speed along the wire towards 
the other lead, where it is also absorbed. The final wire (bottom curve) has a conductance of 
G = 6 Go- Similar thinning by motion of a kink has been reported in TEM experiments.^'^ 

4. UNIVERSALITY OF THE EQUILIBRIUM SHAPE 

Radius functions R{z, t -^ oo) in the interval [0, L/2\ for equilibrium shapes obtained from 
14 simulations of random wires are shown in Fig. 4(a). The conductance of those wires ranges 
from 1 to 200 6*0 with lengths 200 ^ kpL ^ 600. They all have in common a cylindrical 
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Figure 3. Thinning of a wire through nucleation and subsequent motion of a kink. From 
top to bottom, the solid lines show -R(-2, 1960 + r) for part of the wire at times indicated on 
the right of each line (curves for r < 150 shifted vertically for clarity), while the dashed line 
shows an intermediate step of evolution during the nucleation of the kink/antikink pair. 
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Figure 4. (a) Equilibrium shape for 14 different initial random wires, with conductances 
ranging from 1 to 200 Go, and length 200 ^ kpL ^ 600; (b) Same wires rescaled by their 
maximum radius Rmax] (c) Delaunay unduloids of various curvature. 



part with a spherical-like lead. Fig. 4(b) shows that rescaling all lengths by the maximum 
radius Rmax of each wire reveals a regularity of the shape of the lead: It is actually a close 
approximation of a Delaunay unduloid,^*' whose curvature is determined solely by the ratio 
of the radius of the cylindrical part and Rmax- Figure 4(c) shows several unduloids of various 



curvature values, to be compared with the rescaled shapes of Fig. 4(b). The Delaunay unduloid 
of revolution is a surface of constant mean curvature and is an unstable steady state of diffusion 
equation (3) without the shell-effect term. In our case, the deep minima of the electron-shell 
potential. Fig. 1, pin the unduloid at its connection with the cylindrical part, thus stabilizing 
it. In fact, if one switches off the electron-shell potential in the simulations, the equilibrated 
wires break apart, as expected from the Rayleigh instability. 

5. CONCLUSIONS 

We have shown that a dynamical model based on the nanoscale free-electron model^^ and 
including surface self-diffusion of a continuous ionic background accounts for the observed 
formation of nearly perfect cylindrical wires in TEM experiments. The thinning of the wire 
appears to occur in very similar ways in both experiments and simulations: The wire necks 
down near one lead, leaving a kink, or step, that moves towards the other lead at more or less 
constant speed. We show that the shape of the equilibrium wire is universal, and consists of 
a cyhnder connected to unduloid-like leads. 
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